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The ACDM framework offers a remarkably good description of our universe with a very small 
number of free parameters, which can be determined with high accuracy from currently available 
data. However, this does not mean that the associated physical quantities, such as the curvature 
of the universe, have been directly measured. Similarly, general relativity is assumed, but not 
tested. Testing the relevance of general relativity for cosmology at the background level includes 
a verification of the relation between its energy contents and the curvature of space. Using an 
extended Newtonian formulation, we propose an approach where this relation can be tested. Using 
the recent measurements on cosmic microwave background, baryonic acoustic oscillations and the 
supernova Hubble diagram, we show that the prediction of general relativity is well verified in the 
framework of standard ACDM assumptions, i.e. an energy content only composed of matter and 
dark energy, in the form of a cosmological constant or equivalently a vacuum contribution. 

However, the actual equation of state of dark fluids cannot be directly obtained from cosmological 
observations. We found that relaxing the equation of state of dark energy opens a large region of 
possibilities, revealing a new type of degeneracy between the curvature and the total energy content 
of the universe. 

PACS numbers: 


I. INTRODUCTION 

The near-Euclidean nature of the universe was one 
of the main predictions of inflationary cosmology [lj. 
Indeed inflationary models actually predict that the 
curvature radius of the universe should be much larger 
than the Hubble scale do (t 0 being the present age of 
the universe), provided enough e-fold occurred. This 
leads to a tight prediction on the spatial-curvature 
parameter < 10 -4 [|] which can only be bypassed 
at the price of unnatural fine-tuning Q. Within the 
ACDM framework, the curvature of our universe can 
be constrained from the observed cosmic microwave 
background (CMB) fluctuations [3J and current estima¬ 
tions of curvature parameter are achieving impressive 
accuracy. For instance Planck collaboration @ reported 
Ulk = 0.0008 ± 0.0004(95%) when the CMB and some 
other cosmological data are combined. These recent 
results appear as a remarkable evidence in favor of 
inflationary cosmology @|, while the value of the cur¬ 
vature radius of the universe appears as an important 
and fundamental characteristics of the universe. The 
successful standard ACDM picture is leading to the 
existence of dominant and unknown dark components: 
the non-baryonic dark matter and the dark energy, 
which is used as the generic term for the origin of the 


* Electronic address: zolniero@lapp.in2p3.fr 

* Electronic address: alain.blanchard@irap.omp.eu 


accelerated expansion of the universe, as revealed by the 
Hubble diagram of supernovae type la [2, [§]. While there 
is almost no doubt on the reality of the acceleration of 
the expansion of the universe, its physical origin eludes 
us. A large variety of explanations have been proposed: 
scalar field domination known as quintessence, gener¬ 
alized gravity theory beyond general relativity (GR) 
or inhomogeneous cosmological models. The present 
day situation is therefore quite paradoxical: within the 
simplest framework of GR and ACDM picture, cosmo¬ 
logical parameters can be determined with a few percent 
accuracy but on the other hand we have little evidence 
on the validity of this framework. It is certainly desirable 
to test the theory at the foundation of our cosmological 
model Q. The remarkable consistency of the present 
universe with a spatially flat space and its importance 
for the theory of inflation, both encourage for a close 
examination of the actual constraints that can be set 
and to their connections with the energy contents of the 
universe. In the standard analyses aimed at measuring 
cosmological parameters, the approach used does not 
offer a measurement of the curvature of space, as would 
provide a purely geometrical test [lOj. As a crucial 
property of our universe, one would like to measure the 
spatial geometry of our space in a model-independent 
way. The knowledge of the luminosity distance with 
redshift does not allow to determine spatial curvature 
of our 3D space El, and should be complemented by 
the measurement of another observational quantity such 
as the Hubble expansion rate history H(z). This can 
be obtained from longitudinal BAO Ei or from the 
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look-back time of the universe 00 . With the latter 
method, a first estimate of the spatial curvature leads to 
a result consistent with flatness but with uncertainties 
of the order of unity 0| . As GR implies that the 
curvature of the universe is related to its energy content, 
an independent measurement of the latter would allow 
to test if GR holds at cosmological scales. Given the 
remarkable accuracy of GR at describing gravity at 
small scales, it is clearly important to check whether 
GR actually holds on cosmological scales. The relation 
between the curvature and energy contents is precisely 
such a test of GR at the background level. In the 
following, we focus on how to examine the status of this 
relation given the quality of relevant cosmological data. 


II. FROM NEWTONIAN GRAVITY TO 
GENERAL RELATIVITY 

The commonly used metric of a space-time in which 3D 
spatial slices are homogeneous is the Robertson-Walker 
metric: 


The dynamical evolution of the expansion factor a(t) is 
directly related to the energy content of the universe and 
to the underlying theory of gravity. It is possible to derive 
the equation governing a(t) in Newtonian cosmology, by 
writing the mechanical energy of a sphere, the radius of 
which can be taken arbitrarily small (the full justification 
relies on the Birkhoff-Jebsen’s theorem). Because the 
equation does not depend on the radius, GR should lead 
exactly to the same equation [22| . This equation reads: 
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where the pi are the energy densities of the different fluids 
contributing to the gravitational field of the universe and 
K an arbitrary constant. 

The normalized curvature parameter f l k = 
— kc 2 /(Ha ) 2 , where H is the Hubble parameter 
a/a, is introduced to quantify the geometry of the 
universe. In a similar way, we can define a dynamical 
quantity: 


^„ = i-E^ c ( 6 ) 


ds 2 = —c 2 dt 2 + a{t ) 2 - —- + r 2 (d0 2 + sin 2 Odq i> 2 )^ 

( 1 ) 

in which a{t ) is the expansion factor of the universe and 
k is a constant related to the geometry of space. A value 
of k = +1 corresponds to a spherical space, k = 0 to an 
euclidean space and k = — 1 to a hyperbolic space. 

The knowledge of the a{t ) function allows to derive the 
coordinate r: 
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where: 
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and thereby the angular-diameter distance: 

a 0 r 


Da{z ) = a(t)r = 


1 + z 


(4) 


where are the density parameters of the various en¬ 
ergy contents of the universe (f^ = 8nGpi/(3H 2 )). In 
a theory leading to Newtonian dynamic for a(t) (eq. 0, 
the relation would be of the form: 

= a ^k geo ( 7 ) 

with a = K/{kc 2 ) being a constant of the theory. 

When derived in the frame of GR, the constant K is 
not arbitrary anymore but related to the geometric con¬ 
stant k found in the Robertson-Walker metric: 

K = kc 2 (8) 

i.e. a = 1, leading to the well-known Friedman-Lemaitre 
equation which now reads: 

O fc = l-E^ c (9) 


III. A SIGNATURE OF GENERAL 
RELATIVITY IN COSMOLOGY 


or equivalently the luminosity distance (Dl(z) = 
Da(z)(1 + z) 2 ). Eq. 0 makes it clear that the knowl¬ 
edge of Da(z) is not sufficient to infer both the curva¬ 
ture ao and the dynamical evolution of the expansion a(t) 
0- This leads to a well-known degeneracy between the 
curvature and a possible evolving dark energy 00 . 
Breaking these degeneracies within GR is possible by us¬ 
ing different observables, such as for instance growth rate 
measurements 0]. Indeed, the dynamic of perturbations 
may provide information which cannot be obtained from 
the evolution of the background 00 , allowing to dis¬ 
tinguish modified gravity from dark energy models (0. 


Alternative theories of gravity to GR offer a potential 
origin for the acceleration of the expansion, and have 
been thoroughly investigated in recent years. As we have 
seen that independent measurements of fl k on one side 
and of the on the other side allow to test the validity 
of GR. This test also holds for other theories which would 
exhibit a different relation between curvature and matter- 
energy content such as the backreaction models 0] or 
some f(R) models. It is clear that these two quantities 
are physically independent and the prediction of GR is: 

^k geo = ^kdyn ■ 


( 10 ) 
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The Hubble expansion factor can now be written as: 

H{zf = HqE(z) 2 (11) 

= Hq [fl m (l + z ) 3 + Qk dyn (1 + z) 2 + Qde(z)] . 

An estimate of curvature of space has been recently ob¬ 
tained [25l ] from distance measurements, independently 
of assumptions about the evolution of dark energy, but 
still within the context of GR. As we mentioned in the 
beginning, methods exist for measuring the curvature of 
our 3D space independently of the underlying theory of 
gravity by using distances and measurements of the ex¬ 
pansion rate H(z), allowing to test the Copernican prin¬ 
ciple by estimating f2/o(z) (26|. H(z) could be inferred 
from the evolution of redshift with time [23,123, b y dif¬ 
ferentiating the age of the universe (13] or by Baryon 
Acoustic Oscillation (BAO) measurements jl5| . 

The situation is more problematic when one wants to 
estimate from the observations. In GR the energy- 
momentum tensor can be evaluated from the measure¬ 
ment of the local quantities pit) and p{t) and their prop¬ 
erties. However, the dark energy fluids are likely to 
escape to local measurements and their properties can 
probably be constrained only through observations at 
cosmological scales. 

Indeed, a combination of geometrical observations (at 
the background level) cannot provide a way to disentan¬ 
gle the total energy content in distinct components such 
as pressureless matter and dark energy on the other side. 
In the absence of local measurements of the content of 
the universe, densities and pressure, we are left with the 
option to parametrize the observables [23}. 

In the following study, we therefore adopt a 
parametrized description of the dynamical contents of 
the universe and examine whether the prediction of GR 
is verified. 

Dark energy fluids are often characterized by their 
equation of state w(t) relating pressure and density: 
p(t) = w{t)p(t). The simplest dark energy model is cer¬ 
tainly a cosmological constant or equivalently a vacuum 
contribution, corresponding to w(t) = — 1. Therefore, we 
investigate this simple scenario where matter and vacuum 
make up the contents of the universe and then examine 
the case of a dark energy component with a constant w. 


IV. DATA SETS AND METHODS USED 


computes the number of chain steps and the minimum 
X 2 inside it. 

The three standard cosmological observables, SNla, 
large surveys of galaxies and CMB fluctuations are used 
by the program. 

The supernovae set is the JLA sample composed of 740 
SNla resulting from a combination of the first three years 
of SNLS and the three seasons of the SDSS-II completed 
with 14 high z SNla from HST and several low z sam¬ 
ples of SNla [30]. For this work only the diagonal terms 
of the error matrix are used. The expression of the lu¬ 
minosity distance, Dl is modified to introduce the new 
parameters: 


D l 


c{l + z) o 
H 0 



( 12 ) 


The BAO resulting from the couplings between the grav¬ 
itational forces and electromagnetic forces in the pri¬ 
mordial plasma produce a specific pattern in both (the 
fluctuating part of) the radiation field and the matter 
density field. These features provide a powerful tool 
for constraining the cosmological parameters. The con¬ 
straints coming from BAO and CMB fluctuations are 
established with the help of reduced parameters pub¬ 
lished by different collaborations and authors. In prac¬ 
tice we use the X%ao presented in the final publica¬ 
tion of the WMAP collaboration [3l}. This X%ao * s 
built using the values of the quantity r s (zd) / Dv(z mes ) 
or Dy (zm.es ) / r a (zd) at different redshifts z mes published 
by 6dFGS [tI, SDSS-DR7-rec (H, SDSS-DR9-rec[p 
and WiggleZ [35}. The quantity Dy(z m es ), combination 

of radial and transverse distances is defined by the rela- 

l/3 

tion DyiZmes) = [cz(l A Z mes ) D^(Zmes) jH (z mes )] 
at redshift z mes where Da is the angular-diameter dis¬ 
tance between today and redshift 2 and r s (zd), the physi¬ 
cal sound horizon at the end of the drag area (redshift Zd)- 
r s {zd) is given by the relation r s (zd) = f c s (u)/H(u)du 
with c s (u) the sound speed in the primordial plasma. 
The sound speed at redshift 2 is related to the ra¬ 
tio of the baryon density, pb, over photon density p 7 , 
Rb(z) by the relation c a (z) = c/ ^ 3(1 + Rb(z)) where 
R b (z) = 3p b /(4p 7 ). 

The value of r s is computed in our approach using 
the expression [TT] of H(z) with the addition of an extra 
term for the contribution of the relativistic particles at 
the time of radiation matter equality. 


The constraints on the two new parameters of the 
model, flk geo and D k dyn , are obtained by a minimization 
of total likelihood x 2 with the help of a custom Markov 
chain program. We start with eight random points in the 
parameter space and test the conver gen ce of the chains 
with the Gelman-Rubin diagnostics [29[ which is based 
on the comparison of the variance of the parameters in¬ 
side and between the chains. When convergence is estab¬ 
lished, the parameters of each step are stored. For each 
cell in the Qk geo — fi k dyn parameter space, our program 


The statistical properties of the CMB fluctuations 
which are a powerful source of constraints on cosmo¬ 
logical parameters, can be summarized through the fol¬ 
lowing three quantities: the shift parameter, R(z *) = 
v'n m H 2 o(l + z*)Da(z*)/c, the acoustic scale, Ia(z *) = 
7r(l + z*)DA(z*)/r s (z *), and the baryon density flbh 2 
where z* is the redshift of the decoupling between matter 
and radiation. We use the values of these 3 parameters 
and the associated covariance matrix provided by [36} us¬ 
ing the Planck archive data. 
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FIG. 1: 1, 2 and 3 a contours in the ttk geo — ^k dy „ plane 
for models with a dark energy component w = — 1. The line 
corresponds to the prediction of GR. In the Hubble expansion 
factor H(z), the parameter is tightly constrained by the 
Planck shift parameter R{z„). 


The fitting formulas of Hu-S ugiyam a and Eisenstein-Hu 
are used to compute Zd et z* mm- The use of these fit¬ 
ting formulas is justified by the fact that these scales are 
determined by the physics of the matter-baryons-photon 
fluids, well understood in terms of atomic physics and 
linear perturbations [13]. These fitting formulas fully ap¬ 
ply as long as the dark energy fluid contribution remains 
negligible prior to the drag period z > Zd- 


V. DISCUSSION AND CONCLUSIONS 

We have established a simple approach to test GR on 
cosmological scales at the background level. This ap¬ 
proach provides a new type of constraints on the cur¬ 
vature of space without assuming GR (but requiring 
assumptions on dark energy evolution) and separates 
clearly the geometrical and the dynamical curvatures. 
Using the three standard probes, we found interesting 
constraints when assuming a dark energy component in 
the form of a vacuum energy (w = —1): 

J n kgeo = -0.0058 ± 0.0052 (la) 

= -0.043 ±0.034 (la) 

Our estimate on Ofc geo and more importantly its uncer¬ 
tainties, are close to the ones obtained by the Planck 
collaboration for Q k , while the uncertainties on our pa¬ 
rameter £l kdyn are noticeably worse. 

At this level, however, from our analysis we found 
that GR is entirely consistent with existing data with 
little room for any alternative. This offers a new evi¬ 
dence for the remarkable agreement between the stan¬ 
dard GR±ACDM model in view of and modern data in 
Cosmology. It is quite remarkable that introducing a 



FIG. 2: 1, 2 and 3 a contours in the O. kgeo — fl k dyn y for 

models having a dark energy component with an arbitrary 
constant equation of state parameter w. The contours are 
marginalized for —2.5 < w < 0.5. The addition of a new 
degree of freedom, w , widens the contours revealing a degen¬ 
eracy between £lk dyn and f 1k gBO 


new degree of freedom (i.e. one parameter extension of 
the ACDM model) does not change the basic figures on 
the properties of our universe in its simplest version. 
The situation changes drastically when the assumption 
on dark energy is relaxed. For example, allowing for an 
additional dark energy fluid with w = —1/3 in the con¬ 
tents of the universe will produce a degeneracy with the 
actual curvature and will not allow to test directly GR 
at the background level with this approach. 

With an arbitrary constant equation of state parameter 
w we found: 

( fl k eo = -0.025 ±0.03 (la) 

\ V kdyn = -0.45 ±0.6 (la) 

{ w = -0.75±g;i| (la) 

This considerably enlarges the possibilities with a signifi¬ 
cant room for alternative theories to standard GR: there 
is a significant degeneracy between curvature and energy 
contents. This may appear somehow not surprising as 
the determination of the energy contents suffers from 
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a degeneracy problem that has been discussed several 
times in the past. However, the situation is quite 
different from the standard model as the same additional 
freedom on dark energy w does not modify much the 
constraints. For instance from WMAP7, BAO and SN 
one gets w ~ —1.0 ± 0.06 and Clk = —0.006 ± 0.007 
[40j |. On the other hand, the limits on f lk geo that we 
obtained are similar to those obtained in GR based 
models without any assumption on dark energy p5^ and 
similar to those based on mo del-independent methods 
[iTtiffl . While similar, our limits on 17*. clearly 
depends on our assumptions on dark energy: constant 
equation of state w (and the speed of sound equal to 
unity). The conclusion is quite interesting: relaxing the 
assumption on dark energy leads to much wider allowed 
range for flk dyn with an equation of state for dark 
energy which is not well constrained, —1.3 < w < —0.5 
(1 er) and with an energy content which is poorly 


known: —1.105 < Qk dyn < 0.15 (1 a). It is reasonable 
to think that more freedom on dark energy will lead 
to even looser constraints, a question needing further 
investigation . 

It is quite remarkable that the standard picture 
(GR+ACDM) fits so well the present day data while at 
the same time we do not have enough independent data 
yet to clarify the nature of dark energy and therefore the 
appropriateness of the theory of GR. The true universe 
may well be much more complicated. This provides 
further motivation for a deeper investigation of the very 
nature of dark energy. 
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